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Abstract. 

We characterize the soliton solutions of the nonlinear Schrodinger equation on the half Une 
with linearizable bovmdary conditions. Using an extension of the solution to the whole line 
and the corresponding symmetries of the scattering data, we identify the properties of the 
discrete spectrum of the scattering problem. We show that discrete eigenvalues appear in 
quartets as opposed to pairs in the initial value problem, and we obtain explicit relations for 
the norming constants associated to symmetric eigenvalues. The apparent reflection of each 
soliton at the boundary of the spatial domain is due to the presence of a "mirror" soliton, 
with equal amplitude and opposite velocity, located beyond the boundary. We then calculate 
the position shift of the physical solitons as a result of the nonlinear reflection. These results 
provide a nonlinear analogue of the method of images that is used to solve boundary value 
problems in electrostatics. 
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1. Introduction 

One of the hallmarks of integrability of a nonlinear evolution equation (NLEE) is the existence 
of exact A^-soUton solutions. It is well known that each soUton is associated to a discrete 
eigenvalue for the scattering problem for the given NLEE via the inverse scattering transform 
(1ST). This has long been known to be true for initial value problems (I VPs) posed on an 
infinitely extended spatial domain (e.g., see |1, 3 J). Recent developments on the 1ST for 
initial-boundary value problems (IB VPs), however, have shown that the same statement also 
applies for problems posed over a semi-infinite interval [4, 6-8, 10, 16, 17]. The purpose of 
this work is to characterize the soliton solutions of IB VPs for integrable NLEEs. 
For concreteness, we consider the nonhnear Schrodinger (NLS) equation 

iqt + qxx-'2v\qf-q = 0, (1.1) 

where as usual v = -1 and v - 1 denote respectively the focusing and defocusing cases. The 
1ST for (1.1) on -oo < x < oo was formulated in [19], where the soliton solutions in the 
focusing case were also obtained, including the well-known one-soUton solution 

^s(x,f) = Ae'f''^+(^'-^'>'^^lsech[A(x-2W-^)], (1.2) 

where k = (V + iA)/2 is the discrete eigenvalue. The IBVP for (1.1) on < x < oo with 
homogeneous Dirichlet or Neumann boundary conditions (BCs) at the origin was studied in 
[4] using the 1ST on the whole line and an odd or even extension of the potential, respectively. 
The case of homogeneous Robin BCs, 



qMt)-aq{0,t) = 0, 



(1.3) 
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with or e R, was also studied in [5, 10, 18] using a clever extension of the potential to the 
whole line. Recently, a new spectral method was proposed for the solution of IBVPs for 
integrable NLEEs [11-13]. The method applies for generic BCs. In general, the solution 
of the IBVP requires solving a system of coupled nonlinear ordinary differential equations 
involving the spectral parameter in order to eliminate the unknown boundary data. The 
method, however, also identifies a class of linearizable BCs. These Unearizable BCs, which 
for the NLS equation coincide with (1.3), allow one to completely linearize the problem and 
to express the solution of the IBVP as effectively as for the IVP. 

Importantly, in all of the above methods the relation between soUtons and discrete 
eigenvalues that exists in the IVP is preserved in the IBVP, yielding solutions of the form 
(1.2). This leads to an apparent paradox, however, since (1.2) does not satisfy the BCs (1.3). 
A further paradox is that numerical solutions of the IBVP for (1.1) show inequivocally that 
solitons are reflected at the boundary. But the soliton velocity is the real part of the discrete 
eigenvalue, which does not change in time. As we show below, the resolution of these 
paradoxes is that discrete eigenvalues in the IBVP appear in quartets, as opposed to pairs 
in the IVP. This means that, for each soliton in the physical domain (in our case, the positive 
X-axis), a symmetric counterpart exists (i.e., on the negative x-axis), with equal amplitude 
and opposite velocity, whose presence ensures that the whole solution satisfies the BCs. The 
ostensible reflection of the soliton at the boundary of the physical domain (here x = 0) then 
corresponds simply to the interchanging of roles between the "true" and "mirror" solitons. 

It is worth noting that the method to obtain soliton solutions for the IBVP on the half line 
is similar in spirit to the method of images that is used to solve boundary value problems in 
electrostatics [15]. Here, however, unlike the case of electrostatics, the reflection experienced 
by the soUtons comes accompanied by a corresponding position shift, which is a reminder of 
the nonlinear nature of the problem. 

The outline of this work is the following. In section 2 we discuss the 1ST for (1.1) on the 
half Une with linearizable BCs, and we derive the symmetries of the discrete eigenvalues 
for the scattering problem. In section 3 we obtain the precise relations between discrete 
eigenvalues and norming constants. In section 4 we discuss the behavior of the solitons and 
we compute the shift originating from the reflection at the boundary, showing that this shift 
depends on the BCs. Finally, section 5 concludes with some final remarks. 



2. Soliton solutions of the NLS equation on the half line 

Consider the IBVP for the NLS equation on the half line with linearizable BCs; that is, (1.1) 
on < X < oo and f > and with (1.3) given. When a = or a — > oo, the BCs reduce to the 
Dirichlet or the Neiunann BCs: q{Q,t) = and qxiOj) = 0, respectively. It is well-known that 
the NLS equation is the compatibility condition of the matrix Lax pair [19, 10] 

Hx-ik[(T3,fi]=Qn, n, + 2ik^[o-3,fi] = Hn, (2.1) 

where [A, B] = AB - BA is the matrix commutator, and 

H(j:,r.*) = -,QQ<r3-,Q,o-3-2tQ = (_|^"^'52j^ (2-^> 

with r(x,f) - vq*(x,t) and where the asterisk denotes complex conjugation. As usual, we 
assume that qix, 0) is sufficiently regular and decays sufficiently fast as x ^ oo that Jost 
solutions and other relevant quantities are weU-defined. 
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2.1. 1ST for the NLS equation on the whole line 

In sections 2.2 and 2.3 we characterize the solution of the IB VP for the NLS equation on the 
half line using the 1ST for the problem on the whole line and an appropriate extension of the 
potential. Here we therefore briefly introduce the relevant quantities that will be used later. 
We refer the reader to [2, 3, 6, 19] for all details. 

We define the Jost solutions of (2.1) as the simultaneous solutions of both parts of the 
Lax pair (2.1) that reduce to the identity matrix as x — > +oo: for all A: e R, 

H^^\x,t,k) = I + / e*(^-^'>'^3Q(y^;)^(i)(y^;^;t)e-<*(A:-;/)cr3 , (2.3a) 

— oo 

DO 

tP\x,t,k) = I - / e'*(^-^')'^3Q(x',0/^*^^(x',?,A;)e-'^(^-^')'^3 . {2.2,b) 

X 

We denote the columns of the eigenfunctions as iJ-'"\x,t,k) - {f/'"'^\i/'"''^^), m = 1,2. The 
regions of analyticity of the Jost solutions, as effected by the exponentials in (2.3), are [2, 19]: 

^(i,L)^^(2,fi). inifc<0, /z<i'^>,;u(2'^): lm/t>0. (2.4) 

Since detyL/('")(x,f,yt) = 1 for }7i = 1,2, for all e R iJ.^^Hx,t,k) and ^PHx,t,k) are both 
fundamental solutions of (2.1). Hence, for all e R, we have the scattering relation: 

li^\x,t,k) = ii^^\x,t,k)e^''^k{k)&-'^''^ , (2.5) 

where 6{x,t,k) - kx-2k^t. The limit of (2.5) as x — > oo yields, using (2.3), an integral 
representation for the scattering matrix k{k): 

oo 

m = I + / Q-'^^''-^^''^•^^Q{x,t)^^^^\x,t,k)e^'''-^''^'^''^ dx, (2.6) 

— oo 

In turn, (2.6) can be used to establish that the elements a\\{k) and a22(k) of A(k) can be 
analytically continued on hnk < and Imk > 0, respectively. Note that with the above 
definitions, A{k) is independent of time. 

The eigenfunctions and scattering coefficients obey the following symmetry relations: 

H^"''^\x,t,k) = cryfx^"''^\x,t,ky, ti^'"'^\x,t,k) = v(TyH^"''^\x,t,ky, {2.1a) 
for m= 1,2, where 



together with 



1 
V 



aiiik) = a\^{k*) , aim = va\m . i^.lb) 

Equations (2.7) hold for aU values of k for which all terms are well-defined. As a result, one 
can write the coefficients of k{k) as 

Discrete eigenvalues occur when a\i{k) = or aiiik) = [2]. Since in the defocusing 
case V = 1 with vanishing BCs at infinity there are no discrete eigenvalues [19], whenever we 
discuss the discrete spectrum we implicitly set v = -1. Assuming that aiiik)a22ik) VA: g R, 
there exist a finite number of such zeros, since the scattering coefficients are sectionally 
analytic. We also assume such zeros are simple. Let us denote kj for j = 1,...,/ and kj 
for j = I,..., J the zeros of aiiik) and anik), respectively, where Imkj > and Im^^ < 0. 
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The asymptotic behavior of the Jost solutions impUes that discrete eigenvalues occur when 
the decaying eigenfunctions as x — > -oo are proportional to those as oo. That is, 

fi^''^\x,t,kj) = bj^^''^'''*'^n^^''-\x,t,kj) , n^^''^\x,t:kj) = bje-^''^'''''^J^i2^^'''\x,t,~kj) . (2.9) 
One then obtains the following residue relations: 

= Cje^''^''''''jY^''-\x,t,kj), Resf^l =C,e«^'*^V'^ (2.10a) 



Res 

k=ki 



an 



a22 
where 

Cj = bj/d22ikj) , Cj = bj/dnikj) , (2.10fc) 

and the overdot denotes differentiation. As usual, Cj and Cj (or equivalently bj and bj) are 
referred to as the norming constants. The symmetry relations (2.7) imply 

/=/, = bj = -b*, Cj = -C*. (2.11) 

To recover the potential from the scattering data one uses (2.5) to define the matrix 
Riemann-Hilbert problem (RHP) 

M^ix,t,k)-M-ix,t,k) = M^ix,t,k)Rix,t,k), (2.12) 

for all A: e R, where the matrix-valued sectionally meromorphic functions M*(x,?,A:) are 

mx,t,k) = U^'^\xJ,k) /'''')'':''''\ M-(x,?,A:) = (^^^^^^g4^,/z(2^)(x,?,A;)l, (2.13) 
\ a(k) I \ a*(k*) I 



^ix,t,k) = \ n^^M, (2.14) 



a{k) I \ a*(k*) 

the jump matrix is 

v|p(fe)p Q2ie(x,t,k)p(j^) 

and the reflection coefficient is p{k) - b(k)/a(k) for all g R. Since n'''"\x,t,k) -\ + 0(l/k) 
as A: -> 00 for m = 1,2 [2], the RHP (2.12) is solved using standard Cauchy projectors, after 
regularizing by subtracting the pole contributions from the discrete spectrum: 

Mix,t,k) = 1+11 f ^— Res [M+] + ^— Res [M"]) + — ? M+(x,?,fe^) ^'''^'^'^ W . (2.15) 

j=\\k-kjk=kj k-kjk=k'. I 2ni _^ k'-k 

The asymptotic behavior of M'^ix,t,k) then yields the reconstruction formula for the potential 
as qix,t) = -2jUm;fe^ooA:Mj2(^,^,A): 

qix,t) = -2i 2 CjS^'''^''-'*i^fif^(x,t,kj) + - / e^'^^''''*^p{k)nfl(x,t,k)dk. (2.16) 

Hereafter we will often write the norming constants as Cj - Aje.^i^i^^^^i^"!'^'' for j - 1, ... ,7. 

In the reflectionless case [i.e., p{k) = Vfe e R] with v = -1, (2.15) reduces to an algebraic 
system. In particular, if 7 = 1, withAi = (y + /A)/2 and Ci = Ae'*^"'"''^'^^'^^^', one recovers the 
one-soliton solution (1.2). In the general reflectionless case with / > 1, (2.15) and (2.16) yield 
[taking into account (2.10) and (2.13)] the pure multi-soliton solution of the NLS equation as 

where Z = (Zi, . . . ,Zj)^ solves the algebraic system of equations 

(l-i-G)Z = c, (2.17fc) 

with G = (Gjj), c = (ci,...,c/)^, the superscript T denotes matrix transpose and 

J C* e~^''"^*''''^P 

Gjj = Q^'<^Mj)c. 2 ^ , Cj = -2je2'^(^'''*^^C; , (2. 17c) 

p=l {Kj Kp)(Kp Kj') 

for all 7,/ = 1,..., 7. 
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2.2. Dirichlet and Neumann BCs 

When homogeneous Dirichlet or Neumann BCs are given, it was shown in [4] that the IBVP 
for the NLS equation can be solved via the 1ST for the whole line using an odd or even 
extension of the potential. In the case of Dirichlet BCs, one introduces the odd extension of 
Q(x)in (2.1) as 

Q'""(x) = Q(x)0U) - Q(-x)0(-x) (2.18) 

for all X e R, where 0(x) is the Heaviside theta function, defined as 0(x) = 1 for x > and 
0(x) = for X < 0. Then one defines //*^'(x, t,k) through (2.3) with Q(x) replaced by Q'^'"(-^)- 
It is straightforward to see that the following symmetries hold for the Jost solutions and, 
consequently, for the scattering matrix: Vx e R and Vfc e R, 

^i'^^\x,t,k)=^i^^\-x,t,-k), A(-yt) = A-'(yt). (2.19) 

In particular, an(-k) = aiiik). From the above symmetry and (2.7), it then follows that 

a*i-k*)=aik) and b{k) = -b{-k). (2.20) 

Recall that / is the number of discrete eigenvalues in either the upper-half or lower-half 
complex ^-plane. The above relations imply that the discrete spectrum has the following 
properties: (i) J is even; (ii) the discrete eigenvalues appear in quartets, namely 

{+kn,±k*X^,, (2.21) 

so that the number of discrete eigenvalues in each quadrant plane is N = J/2. Moreover, 
(iii) the above synmietry in the discrete eigenvalues also induces a relation between the 
corresponding norming constants, as we show next. Throughout this work, we label the 
discrete eigenvalue symmetric to k„ as 

k„' = -k*„. (2.22) 

We also take Rek„ for n = 1,...,A^ to avoid singular cases. Using (2.7) and the 
symmetries (2.19), we have, for « = l,...,N, 

b„bl,=-l, CnCl,^\la\kn). (2.23) 

The relations between the norming constants take on a particularly simple form in the 
reflectionless case. In particular, for reflectionless potentials and in the case N - 1 the 
trace formula yields simply 0(^2) = W[-A(A + iV)], where ki = (V + iA)/2 (see section 3 for 
details). Writing C„ = A„e'*"^"+'**'""^'^^^) forn = 1, . . .,2N, we then obtain the following exphcit 
relations: 

^i+^2 = ^log(l + ^), cp2-cpi=2aigiA + iV). (2.24) 

Next, consider the IBVP for (1.1) with Neumann BCs. In this case one can use the even 
extension of the potential: 

Q^'"(x) = Q(x)0(x) -I- Q(-x)0(-x) (2.25) 
Vx e R. As before, one obtains symmetry relations for the eigenfunctions and scattering data: 

fi'-^\x,t,k) = cr3i/^\-x,t,-k)cr3 , k{-k) = cr3A-\k)cr3 . (2.26) 
For the scattering coeflficients we therefore have 

a*i-k*) = a{k) , b{k) = b{-k) . (2.27) 
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Thus, (2.22) still applies. Moreover, one can again obtain relations between the discrete 
eigenvalues and the norming constants: 

bnbl,^\, C„Cl,^-l/d\k„), (2.28) 

forn = l,...,N. In particular, for pure sohton solutions and A? = 1 it is (cf. section 3): 

^i+^2 = ^log(l + ^), ^2-^1 =2arg(A + jy) + ;r, (2.29) 
where C„ = A„e^"f"+'((P"+'^''2) ^^foj-^,, 

2.3. Robin BCs 

It was shown in [10] that, even in the case of Robin BCs (1.3), one can still reduce the IB VP 
to the solution of an IVP on the whole line. Motivated by the linear problem (see Appendix), 
one introduces the extension of Q{x) as 

Q^^\x,t,k) = Q{x,i)®{x) + ?{k)Q{-x,t)®{-x) x e R, (2.30) 

where F(A:) = diag(/(A:),/(-A;)), with 

f{k) = (2k - ia) I {2k + ia). (2.31) 

(The factors of 2 in (2.31) are chosen for consideration of the Unear Umit, discussed later.) 
Since F(^) - I when o; = 0, and fik) ^ -I as o- ^ oo, (2.30) reduces to the odd/even extensions 
of the potential in the case of the IB VP with Dirichlet/Neumann BCs, (2.18) and (2.25) 
respectively. More generally, the extended potential Q^^^{x,t,k) satisfies the symmetries 

Q^''\-x,t,k) = F{-k)a'''\x,t,k), (2.32a) 
'L{k)f{-k)Q''''\x,t,k) = -Q^''\x,t,k)I.(k), {2.32b) 
where S(A;) = diag(/(A:), -1). 

Let i/^\x,t,k) for /' = 1,2 be the Jost solutions defined via (2.3) with Q(x,f) replaced by 
Q^^^ix,t,k), as before. Additional care is necessary for the case of Robin BCs compared to 
Dirichlet and Neumann BCs, because, for all a^O: (i) Q^^^ix,t,k) is discontinuous at x = 0; 
(ii) for all x < 0, Q'""(x,f,fc) has simple poles at k - ±ia/2 [since F(k) does]. Similar issues 
arise in the linear limit (cf. Appendix). [There, an extra term must be added to the f-part of the 
Lax pair (A.2) to restore its compatibility at x = 0. That is not necessary here, since the Lax 
pair (2.1) is multiplicative, while that of the linear Schrodinger equation is additive.] Note 
also that, when r(x,t) = vq*{x,t), the compatibility of the extended Lax pair (2.1) with Q(x,/) 
replaced by Q^^\x,t,k) implies that q^^\x,t,k) solves the NLS equation for all x 0, since 
f(k)f(-k) - 1. As in the linear limit, even though Q^^^(x,t,k) is discontinuous at x = 0, the 
BCs (1.3) are satisfied from the left and from the right. 
From the synametries (2.32) one obtains 

ix'-^\x, t,k) = I.ik)fP\-x, t, -k)I,i-k) , (2.33a) 
Ai-k) = 'L{-k)k-\k)'L{k) . {2.33b) 

In particular, 

a{k) = a*i-k*), b(k) = f(k)b(-k) . (2.34) 

Equations (2.33) and (2.34) generalize the relations found for Dirichlet and Neumann BCs 
[namely (2.19), (2.26) and (2.20), (2.27), respectively.] 

Regarding the analyticity of the eigenfunctions, note first that i/^\x,t,k) is the same as 
that of IVP for all x > 0, and therefore enjoys the same properties there. Using (2.33a) and 
definition of one can obtain the regions of analyticity of ^^^\x,t,k) for x > 0: 
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• ju*^'^' is analytic for ImA: < 0, except for a simple pole at A: = iajl when or < 0; 

• /i*'"''^ is analytic for \mk > 0, except for a simple pole at A: = -jar/2 when or < 0. 

[The analyticity properties for x < can be obtained using the symmetries (2.32), (2.33) and 
the integral equations (2.3). These properties are not necessary for our purposes, however, 
since we are only interested in reconstructing the extended potential for x > 0.] The integral 
representation of the scattering matrix obtained from (2.5) in the hmit x — > oo and the 
symmetry (2.32) imply 

oo oo 

a{k) = 1 - f{k) / il^li-x, t, -k)r{x, dx - / )uf/(x, t, -k)r{x, t) dx . (2.35) 



Thus, a{k) is analytic for ImA > 0, except for a simple pole at A = -iajl when a <0. Thus 
lJ-^J^\x,t,k)/a(k) in (2.12) has a removable singularity at A = -/a/2. By symmetry, the 
same apphes for ^^^'^\x,t,k) I a* {k*) in the lower-half plane. Therefore, the RHP defined 
in section 2. 1 also applies to the case of Robin BCs (1 .3), and for all x > 0, the potential q{x, t) 
is also reconstructed in the same way. 

Equation (2.34) implies that the symmetry property (2.22) of the discrete spectrum 
applies in the case of Robin BCs as well. Moreover, similar symmetries as before exist for the 
nomaing constants (again, cf. section 3): 

for n = l,...,N. In particular, for pure soliton solutions with A' = 1, the following relations 
exist between the norming constants associated to symmetric eigenvalues: 

1 . /. 1 . ry2 + (A-a)2- 



bnbl, = f(kn) , CnC*, = - , (2.36) 



^l+^2 = ^l0g(l + ^)+-l0g 



(p2-ipi =2arg(A-l-jy)-arg 



V + i(A-a) 



(2.37a) 



+ 71. (237 b) 



V + i(A + a)\ 

Equations (2.36) and (2.37) reduce respectively to (2.23) and (2.24) as a — > oo and to (2.28) 
and (2.29) as a ^ 0. In section 4 we discuss what (2.22) and (2.36) imply about the physical 
behavior of the soUtons in the IBVR 



3. Relations between discrete eigenvalues and norming constants 

We now derive the relations between the discrete eigenvalues and the norming constants 
we presented in section 2, together with their generalization for multi-soliton solutions and 
solutions with a nonzero reflection coefficient. We discuss explicitly the case of Robin BCs 
[namely (2.36) and (2.37)]. The corresponding relations in the case of Dirichlet and Neumann 
BCs [namely, (2.23), (2.24), (2.28) and (2.29)] follow from (2.36) and (2.37) by taking the 
limits a ^ oo and a 0, respectively. 

Recall first that the discrete eigenvalues and norming constants obey the synnmetries 
(2.11) and (2.21). From (2.7) and (2.33a) we have, for all n = 1, ... .AT, 

H^^'^\x,t,kn>) = b„,e^'''^''''*"'^n^^''^\x,t,k„'), (3.1a) 

as well as 

/|(i'^>(x,f,A„0 = -fo-ie2«(-^''-^")2(A„0//^''^\-x,f,A„), (3.1^) 
where A„' = -A* as before. Also, using (2.33a) we can write 

fx'-^'^\x,t,k„') = f{knW\kn)tP'^\-xM . (3.2) 



Solitons, EVPs and a nonlinear method of images 



8 



Inserting (3.2) into (3.1) leads to 

^ikn')fl^^'^\-X,t,k„) = -bnbn'f{knW\kn)tl^'^\-X,t,kn) ■ (3.3) 

Since f(k„) = 1 1 fik„>) and f{k„>) = f*(k„), we then have fo„fo„' = -f*{k„) for all n = 1, . . . .A'', 
which in turn, using (2.11), yields the first of (2.36). 

Now recall the definition of the norming constants (2.10) and the symmetry (2.1 1). Using 
the symmetries of the scattering coefficients il.lb) and (2.34), we obtain cinikn) = —iijii—kn) 
for all n = 1, . . .,2N. Hence C„ = bl/d22ikn'), and then, noting that d{k„') = -d*(kn), we obtain 
the second of (2.36). Note also that, when Cj = Aj&^iij^'^fi^''!^\ the second of (2.36) impUes 

^„+^n' = (log|/(A:„)|-21og|fl(A:„)|-21ogA„)/A„, (3.4fl) 

Vn - ^n' = arg[/(/:„)] - 2 arg[a(fc„)l + n . (3.4fo) 

We now derive (2.37) and its generalization to arbitrary solutions of the IB VP. It is well 
known that the analytic scattering coefficients obey trace formulae. Explicitly, for the NLS 
equation, a{k) = aiiik) is given by [2]: 

lo^aik) = 2 log (— ) + - / d^ , (3.5) 

for all ImA: > 0. Using (3.5) in (3.4) then yields half of the norming constants in terms of the 
other half. In particular, for reflectionless solutions the integral in (3.5) vanishes, and (3.5) 
yields simply 

d{kj) = n ' ikj -km) In ikj - kl,) (3.6) 

m=\ ' m=l 

for all j = I,..., J, where the prime indicates that the term with m - j is omitted from the 
product. Using the symmetry of the discrete eigenvalues, (3.6) becomes 

^^j^ ^ _ kn+k^ ^ , (kn — km){kn + A^) 
Ikfiikfi ~ kfi) fn=\ (.kn kfti)(kn — k^) 

V„ ^ , [ V„ - V,„ + i(A„ - A,„)] [V„ + V„, + i(A„ - A,„)] ^^^^ 



iAniVn + iA„) „=i [Vn + Vm + iXA„ + A„)][y„ - + i(A„ + Am)] ' 
for all n = l,...,N. One can now substitute (3.7) into (3.4) to obtain the generaUzation of 
(2.37) as: 

1 . f. Al\ 1 . rV^ + (A„-a)2^ 



Vi + (A„+ay 



± l(Vn - Vmf + (An -A^f] [(V„ + Vmf + {An - Amf] 

'An „=1 [(Vn + y„)2 + (An + A„)2][(y„ - Vj^ + (An + A^)2] ' 



- = -2 arg {An + i y„) + arg 



Vn + i{An-a) 



Vn + i{A„+a) 



-n 



^ ^ , [ V„ - + i{An - A„)] [y„ + y„ + i{An - Am)] _ , , 

-2 y arg . (3.8o) 

„tl ^[Vn + Vm + i(An+Am)Wn-Vm + i(An+Am)'\ 

In the simplest case A' = 1 there is no product in (3.7), and the relations between the norming 
constants become especially simple: 

CiCl = {a-A + i V){A - i V)^A^I[{a + A-iV)V^] 
[where ki = {V + iA)l2 as before], which, when inserted in (2.36), yields (2.37). 

Note from (3.4) and (3.7) that, when N > 2, the norming constants of a given mirror 
soliton are affected by the presence of all the physical soUtons. This is due to the fact that 
every soliton interaction produces a position and phase shift, as we discuss next. 
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4. Soliton behavior 

We now discuss the behavior of the soliton solutions of the NLS equation on the half hne. 
We refer to the soUtons located to the right of the boundary (i.e., at x > 0) as the physical 
solitons, and to the counterparts of the physical solitons to the left of the boundary as the 
mirror sohtons, since they can be considered as a reflected image of the physical sohtons, as 
we will see. Equation (2.21) obviously implies that the number of physical solitons equals 
that of mirror solitons, and, with the above notations, this number is A'^ = J 12. 

4.1. Soliton reflection 

We first discuss the case A' = 1 for simplicity. Solving the algebraic system (2.17) one obtains 
the two-soUton solution of the NLS equation, and the solution of the IBVP is then obtained 
by choosing the norming constants of the mirror solitons as explained earlier. 

Let ki = (y + J A)/2 be the discrete eigenvalue of the scattering problem corresponding to 
the physical soUton. Recall from (1.2) that the real part and imaginary part of the discrete 
eigenvalue determine respectively the velocity and the amplitude of the soliton. By the 
symmetry (2.22) of the discrete spectrum, we know that k2 = -k\ - {-V + iA)l2 is also 
an eigenvalue, corresponding to the mirror soUton. Hence, the mirror soUton has the same 
amplitude as, and opposite velocity to, the physical soliton. 

Figure 1 shows (left) the soliton reflection at the boundary in the case of Dirichlet BCs, 
together with a contour plot (right) that includes the mirror soliton (dashed lines). The above 
results imply that the soliton reflection at x = is simply a particular case of an elastic two- 
soliton interaction of the NLS equation, in which the norming constant of the mirror soliton 
is chosen so as to make the whole solution zero at the origin. As usual in a soUton interaction, 
the solitons re-emerge intact after the collision, except that in our case the roles of physical 
and mirror soliton are now swapped. A similar scenario occurs in the case of Neumann and 
Robin BCs, as shown respectively in Figs. 2 and 3, except that the norming constant of the 
mirror soliton in each case is such that the appropriate BCs are satisfied. 

Note that the symmetry of the discrete spectrum and the relations between eigenvalues 
and norming constants apply independently of whether the physical soUton has a positive or 
negative velocity. In other words, they apply whether the discrete eigenvalue associated to 
the physical soliton is located in the first or second quadrant of the complex A:-plane [recall 
V„ = 2Rekn]. Of course, if the physical soliton has a positive velocity, no soUton reflection 
occurs for f > 0, and the solution is exponentially small at the origin for all f > 0. Nonetheless, 
a mirror soliton is still needed to satisfy the BCs at the origin, as shown in Fig. 4. 

The location of the mirror soliton is the same for Dirichlet and Neumann BCs, and in 
those two cases only the phase diflFerence between the solitons contributes to satisfying the 
BCs [cf. (2.24) and (2.29)]. The same is not true, however, for the more general BCs (1.3) 
witiiQr^O[cf. (2.37)]. 

Finally, figure 5 displays the reflection of two physical solitons in the case of Dirichlet 
BCs, demonstrating that our results are not limited to the case N = I. Similarly, figure 6 
displays the reflection of a two-soliton bound state, also in the case of Dirichlet BCs. As 
before, when > 1 one first solves the algebraic system (2.17) to obtain the 2Ar-soliton 
solution of the NLS equation. The solution of the IBVP is then obtained by appropriately 
choosing the norming constants of the mirror solitons. The generalizations of (2.37) to obtain 
the norming constants for solutions with N >2 and for solutions with non-zero reflection 
coefficients are described in section 3. 
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Figure 1. Soliton reflection at the boundai'y in the case of Dirichlet BCs, with A = 2,V = -2, 
f = 5, and ip = 0. Left: three-dimensional (3D) plot of |^(.y, /)|. Right: contour plot showing the 
mirror soliton (dashed) to the left of the boundary. 




Figure 2. Soliton reflection at the boundary in the case of Neumann BCs, with A = 2, V = -2, 
f = 5, and ip = 0. Left: 3D plot. Right: contour plot showing the mirror soliton. 




Figure 3. Soliton reflection at the boundary in the case of Robin BCs with o; = 3, and with 
A = 2, V = -2, f = 5, and ip = 0. Left: 3D plot. Right: contour plot showing the mirror soliton. 
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-5 



Figure 4. Reflection at the boundary in the case of Robin BCs and positive velocity: a = 3, 
A = 2, V=l,^ = 3, and if = 0. Left: 3D plot. Right: contour plot showing the mirror soliton. 
In this case the solution with Dirichlet and Neumann BCs is visually undistinguishable from 
the above. 




Figure 5. Reflection at the boundary of two physical solitons: Ai = 2, At = 3/2, Vi = -3, 
V2 = -1,^1 = 12, ^2 = 8, and if[=(p2= 0. Left: 3D plot. Right: contour plot. 




Figure 6. Reflection at the boundary of two physical solitons: A[ = 1, A2 = 3, Vi = V2 = -1, 
^1=^2 = 10, and ifii =ifi2 = 0. Left: 3D plot. Right: contour plot. 
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It is convenient to label the discrete eigenvalues such that Refci < ReA:2 < • • • < ReA:2iv- 
With this convention, as f — > -oo all physical solitons correspond to the discrete eigenvalues 
in the second quadrant of the complex A:-plane (i.e., those with Re^„ < 0). Then, as each 
soliton is reflected in succession, the corresponding discrete eigenvalue associated with the 
physical soliton switches role with its symmetric conuterpart, until, as f — > oo, all physical 
solitons correspond to the A' discrete eigenvalues in the first quadrant (i.e., with Re/:„ > 0). 
Consequently, as t ^ -oo the discrete eigenvalues associated with the physical solitons as 
ki,...,kN, and the corresponding mirror solitons are given respectively by k2N, fcv+ 1 • That 
is, n' = 2N-n + 1 for all n = 1,...,A^. Conversely, as f ^ oo the eigenvalues associated 
with the physical solitons are k^+i,. . .,k2N, and the corresponding mirror solitons are given 
respectively by kiq,...,k\. That is, n' = 2Ar-n+ 1 for all n = N +l,...,2N. 



4.2. Reflection-induced shift 

Recall that the soliton reflection at the boundary is effectively the interaction between the 
physical soliton and its mirror image. Since any soUton interaction results in a position shift, 
it follows that the soliton reflection at the boundary also produces such a shift. As we show 
next, however, a second contribution also exists to the total reflection-induced shift. In fact, 
depending on the soliton parameters, this second contribution can even make the total shift 
zero as if the whole process were purely linear 

It is well-known that, as f ^ oo, a multi-soliton solution becomes asymptotically a hnear 
superposition of one-soUton solutions [2, 19]. That is, 

q{x,t) ~ q%x,t) (4.1) 
7=1 

as f — > ±oo, where qj(x,t) is of the form (1.2), but with the soliton parameters A, V,^,^ replaced 
by Aj,Vj,^J,(p'^ for j = l,...,J. The fact that ^ and (pj do not coincide with each other is 
the manifestation of the interaction-induced shift. Define as usual the interaction-induced 
position shift of the j-th soliton as - -^J, and label the discrete eigenvalues so that 
Vi <V2 < ■■■ <Vj. Without repeating the calculations [2, 19], we quote the relevant results: 
for all 7 = 1 , . . . , 7 it is 

J J 
^j^j = ^j^] ~ log 1^(^7)1 ~ li ' log \amikj)\ = Aj^j - log \d(kj)\ + i; ' a^j log la^ikj)] , 

m=l m=l 

(4.2) 

where ajik) = (k- kj)/(k- k*) is the transmission coefficient for a one-soliton solution [cf. 
(3.5)], amj = -1 for m = 1,. . - 1 and (Tmj = 1 for m- j+l,...,J, and as before the sum is 
taken over aUmi^ j. Comparing the asymptotic results as ? ± 00, we then obtain the position 
shift for multi-soliton solutions of the NLS equation as 



2 

^^7 = li'o^mjlog 
^; m=l 



k-k* 



(4.3) 



Note that the position shifts are pairwise additive. That is, - Y/m=i ^^j,m^ where 5f ;>„ is the 
position shift of the 7-th soliton arising from its interaction with the m-th soliton. In particular, 
if ^„ = (y„ -H iA„)/2 with V„ < is the eigenvalue associated to the physical soliton at ? = 0, and 
k„' = -k* is the mirror eigenvalue, (4.3) yields the position shift resulting from the interaction 
of the physical sohton with its mirror: 

SUn' = - l0g|l + ^ j = -6^n'^ . (4.4) 
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^2 V, 



Figure 7. Left: schematic diagram of the reflection-induced sliift A^. The displacement from 
the origin of the interaction center has been exaggerated for illustration purposes. Right: a 
soliton reflection with a large, boundary-induced shift: a = -2, A = 2,V = -\,^ = 5 and ip = 0. 



Since the physical and mirror soliton interchange roles in the reflection, however, the 
total reflection-induced shift is not simply due to the soliton interaction. More precisely, define 
the reflection-induced position shift as the displacement of the soliton center from where 
it would be had it followed a purely piecewise linear path [cf. Fig. 7]. A simple calculation 
shows that such a position shift is given by 

^^n=^n+Q-6^n, (4.5) 

where as before 6^„ the interaction-induced shift, obtained from (4.3), and the index «' 
labels the mirror soliton of the n-th soliton. As is evident from (4.2), the asymptotic soliton 
parameters and i^j" do not coincide with the constants and (fj appearing in the exact 
7-soliton solution (2.17). In particular, using the symmetry of the discrete eigenvalues, after 
some tedious but straightforward algebra (4.2) yields simply 



=^„+^„' +21og|fl(fe„)| = log|/(fe„)|-21ogA„ = — log 



Vl + {A„-a) 



2> 1/2 



(4.6) 

where (3.4) was used. The reflection-induced shift A^ is then obtained inserting (4.4) and 
(4.6) into (4.5). In particular, for = 1, it is: 



A^ = 



1 



(4.7) 



,y2 + (A + a)2y 

Importantly, the reflection-induced shift depends on the BCs. Indeed, such a dependence 
is evident in Figs. 1, 2, 3 and 7. It should also be clear than A^ can be either positive or 
negative depending on the soliton parameters and BCs. Note also that, if the physical soliton 
has a positive velocity (i.e., if V > 0), the reflection at the boundary and the corresponding 
position shift obviously occur at f < 0. 

The above results are easily generalized to the case N >2. Indeed, (4.6) holds for all A^. 
Also, taking n < N (i.e., V„ < 0), using the symmetry of the eigenvalues, (4.2) yields 



A„ 



1 + - 



y2 

' n 



1 

An 



N 
m=l 



log 



(Vn-V„y + (A„-AJ^ 



(V„-y„,)2 + (A„-fA„)2 



1 

An 



N 

li 'o-n 
m=l 



,"log 



(V„ + Vm)' + (A„-A„)2 



(V„ + Vm)2 + (A„+A„)2 



(4.8) 



One can then combine (4.6) and (4.8) into (4.5) to obtain the generalization of (4.7). 
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Note that, when A'^ > 2, the interaction shifts 6^n compare the soliton positions before 
and after all interactions have occurred. Similarly, the constants and if^ apply for times 
respectively before and after all soliton interactions (and thus reflections) have occurred. So 
(4.5) applies in that limit. But equations (2.36) allow one to obtain the A^-soUton solution of 
the IB VP given the A' discrete eigenvalues and norming constants associated with the physical 
solitons at any value of t, whichever these may be, and independently of how many reflections 
may already have occurred. 

5. Discussion 

Since many physical situations naturally give rise to IBVPs for the NLS equation, we expect 
that our characterization of the soliton solutions of the IBVP will have a broad range of 
applicability. On the technical side we note that the approach we used for the IBVP — namely, 
extension of the potential and use of the 1ST for the IVP [4, 10] — is fundamentally difl'erent 
from the new method for IBVPs presented in [1 1], which is based on the simultaneous spectral 
analysis of both parts of the Lax pair. We also note that the symmetries of the scattering 
coefiicients had been derived in [4] for Dirichlet and Neumann BCs and in [10] for Robin BCs 
(see also [5, 18]). None of those works, however, discussed the symmetries of the discrete 
spectrum, norming constants and the corresponding implications on the soliton behavior. 

We should also emphasize that the symmetries of the scattering data and the discrete 
spectrum only apply to linearizable BCs. In fact, a trivial counterexample in the case of other 
BCs is given by (1.2), which provides a perfectly valid solution to the IBVP for the NLS 
equation on the half line with the non-homogeneous Dirichlet BC ^(0, t) = qsiO, t). 

Finally, we note that, similarly to the traditional method of images, one could in principle 
assume that the symmetry (2.21) of the discrete eigenvalues holds, require that the solution 
satisfies the given BCs and then obtain the relation between the norming constants by solving 
a system of algebraic equations. We did not do so here, however. On one hand, the resulting 
equations are transcendental, and finding their solution without any a priori knowledge is 
highly nontrivial, which would make the method impractical. Moreover, exploiting the 
symmetry of the scattering coefiicients enabled us to prove that all solutions of the IBVP 
possess the syrmnetry (2.21) and to find the relations between norming constants, regardless 
of the number of physical soUtons and of whether the reflection coefficient is zero. 

The present results open up a number of interesting questions. An obvious one is what 
happens for other integrable NLEEs. We expect that similar results will apply to equations 
such as sine-Gordon equation [9] and the Ablowitz-Ladik lattice [6]. Less clear, however, is 
what happens for the Korteweg-deVries equation, since in this case all solitons travel to the 
right (or to the left depending on the sign of the nonlinear term), and no obvious symmetry 
exists. It might be that the nonUnear method of images only appUes as presented to equations 
that admit some reflection symmetry, or such that the linear limit can be solved with Fourier 
methods [14]. (Note that even the traditional method of images has similar restrictions.) 

Another nontrivial question is what happens on finite domains. Numerical simulations 
show that in this case the solitons experience an infinite number of reflections at each 
boundary. After two reflections, however, any sohton recovers its original velocity. The initial 
soliton and its twice-reflected copy must thus be associated to the same discrete eigenvalue. 
Moreover, since the soliton experiences an infinite number of refiections, the extension of 
the solution to the whole line would require an infinite number of mirror copies. Finally, 
such a solution would not possess decaying or even constant BCs as x — > ±oo. Therefore, its 
characterization seems to be outside the current capabilities of IST-based methods. 
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Appendix: IB VP for the linear Schrodinger equation with Robin BCs 

For comparison purposes, here we solve the IB VP for the Unear Schrodinger (LS) equation 
iqt + qxx = ^ (A.l) 

for X > 0, with the homogeneous Robin BCs (1.3) at the origin, using the 1ST for the IVP 
and an extension of the potential to the negative real x-axis. (See [6, 11] for an alternative 
method.) 

Equation (A.l) is the compatibiUty condition of the scalar Lax pair [6] 

(f>x-ik(p = q, (pt + ik^(p = iqx-kq. (A.2) 
When (A.l) is posed on -oo < x < oo, the lost solutions of (A.2) are 

X oo 

(p^(x,t,k)^ J e*(^-^''^(x',f)dx', <^-(x,f,yt) = -/e'*<^-*'^^(x',Odx'. 

^ (A.3) 

It is easily shown that <p'''(x,t,k) can be analytically continued on 1mA: ^ 0, and <p^(x,t,k) = 
0{l/k) as k CO. Moreover, on 1mA: = one has the following jump condition: 

(P^(x,t,k) - (p-(x,t,k) ^ e"'''q{k,t) , (A.4) 

with q(k, f) given by the first of (A.6) below. Equation (A.4) defines a scalar RHP that is 
solved with standard Cauchy projectors. Moreover, the jump data satisfies 

q(k,t)^e-'''''q(k,0) (A.5) 
and the solution of the RHP leads to the direct and inverse Fourier transform pair: 

oo 1 oo 

qik,t) = / e-'*^^(x,Odx, qix,t) = ^ / e'''''qik,t)6k. (A.6) 

— oo '* — oo 

Now consider the IB VP for (A. 1) on x > with BCs (1 .3). It is trivial to see that, if q{x, t) 
is given by the second of (A.6) for x > 0, the BCs (1.3) are satisfied provided that 

. q{k,t)=fo{k)q{-k,t), (A.7) 

where 

fo{k) = (k-ia)l{k+ia) . (A.8) 

This motivates the following extension of q{x,t): 

q''^\x, t, k) = q{x, t)®{x) + fo{k)q{-x, i)@{-x) , (A.9) 

for all X e R. The corresponding Fourier transform (f^^{k, f) [defined by the first of (A.6) with 
q{x,t) replaced by ^""'(x,?)] is related to the the one-sided Fomer transform qoik,t) 

oo 

§o(A:,0 = /e"%(x,Odx, (A.IO) 



hy q^^\k,t) = qo(k,t)+ fo(k)qo{-k,f). It is then easily verified that ^^''^ (A:,?) satisfies (A.7) [note 
fo{k)fo{-k) = 1]. We then consider an extended Lax pair obtained by replacing q{x,t) with 
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q^^\x,t,k) in (A.2). Note however that q^^\x,t,k) is in general discontinuous at x = 0, and in 
a distributional sense it solves 

iq^"" + qfx = Uo{k) + 1] t)6{x) - Uo{k) - 1] ^(0, t)6'ix) , 

where 6(x) is the Dirac delta. Thus the Lax pair (A.2) with q(x,t) replaced by the extended 
potential is not compatible at x = 0. The compatibility is restored if the f-part of (A.2) is 
replaced by 

4>, + ik^4> = iqf^-kq^^'^ + iifoikJ-lMOjySix). (A.ll) 

One now defines </>'^ix,t,k) as the Jost solutions of this new Lax pair. Since the x-part is 
formally the same as before, they are simply given by (A. 3) with q(x, f) replaced by q^^^{x, t,k). 
For X > 0, (p~ix,t,k) is analytic on Imk < 0. Also, for x > 0, ^^{x,t,k) is analytic on Imk > 
when Rear > and is meromorphic on ImA; > with a simple pole at A: = -ia when Rear < 0. 
[When Rea > 0, the pole at = -ia is in the lower-half plane.] The jump condition (A. 4) still 
holds with q{k,0) replaced by (f^^{k,0). Moreover, using (A.ll) with q{x,t) — > as x — » oo, 
we obtain that the scattering data q^^^{k,t) still evolves according to (A.5). After subtracting 
the pole contribution, one solves the RHP (A.4) using Cauchy projectors, obtaining: 

<l>{x,Uk) = -i, / e'(*'-*"" ^^^^^ dk' + -i- y;. Res [^-] , 
2ni k'-k k + ia " k=-ia^ ^ 

where = for Re a > 0, = 1 /2 for Re a = 0, and v'„ = 1 for Re a < 0. Note that 
Resk=-ia[4>^] = -2iae"^'^'" 'qiia,0). Substituting the above into (A.2) we have, for all x > 0, 

q(x, t)^ — [ e'(*^-^'')^'""(A, 0)dk- 2v'„a&'"'^'''^'qo{ia, 0) , (A. 1 2) 

— oo 

Note that, even though q^^^{x,t,k) is discontinuous at x = 0, the BCs (1.3) are satisfied from 
the left and from the right. 

One can compare the above solution to the linear limit of the solution of the IB VP for 
the NLS equation. To do so, let Q^'"(x,r,A;) = Q^'"(x,?,A;;e) = 0(e) in (2.1). As e ^ 0, it is: 

OO 

a{k) = l + 0{e^), b{k) = / e^''"'q^^\x,t)Ax + 0{e^), (A.13a) 

— OO 

X OO 

n'-l^(x,t,k) = J e^"'^''-'''\ix',t)dx'+m J e^'''^''^''\(x',t)dx' + Ois^). 

(AA3b) 

Since //'^'•'^'(x, t,k) has a pole atk = -ia/2when a < Obuta(A:) = 1 to leading order, the solution 
of the RHP (2.12) acquires an additional term generated by the residue of f/^''^\x,t,k)/a(k). 
Taking the pole contribution into account, and performing the change of variable, 2k —> k', 
one then obtains, to 0(e), (A. 12). That is, in the Unear limit, the solution of the IBVP for the 
NLS equation coincides with solution of the IBVP for the hnear Schrodinger equation. 
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